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Abstract 

We study the Gross-Neveu model in 2+1 dimensions with a baryon chemical poten- 
tial \x using both analytical and numerical methods. For /i greater than a critical 
value the model is chirally symmetric and has a Fermi surface with kp — /i. We have 
calculated the particle interaction in medium due to scalar meson exchange to leading 
order in Nj 1 , where Nj is the number of flavors, in the hard dense loop approach. 
The result has been used to calculate the relation between [i and the Fermi momen- 
tum and velocity in the resulting Fermi liquid to 0(NJ 1 ). Simulation results from a 
32 2 x 48 lattice for fermion and meson dispersion relations and meson wavefunctions 
are then presented, showing qualitative and in some cases quantitative agreement 
with analytic predictions. In particular, the simulations show clear evidence for the 
in-medium modification of the scalar propagator, oscillatory behaviour in the wave- 
function consistent with a sharp Fermi surface, and tentative evidence for a massless 
pole in the vector meson channel resembling zero sound. 

PACS: ll.10.Kk, 11.15.Ha, 11.15.Pg, 21.65.+f, 67.90+z 
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1 Introduction 



Interest in QCD at low temperature with a non-zero background baryon density is 
currently at a high level. Understanding the ground state of any strongly interacting 
many-body system requires non-perturbative techniques, and lattice QCD is usually 
considered one of the most systematic and successful. Unfortunately, attempts to 
develop Monte Carlo simulation techniques for systems with a non-vanishing baryon 
chemical potential \x have so far met with little success [I]; understandably, therefore, 
there has been comparatively little parallel development in relating correlation func- 
tions measured using standard lattice field theory techniques to the phenomenology 
of dense systems. Exceptions are recent attempts to measure diquark condensates 
(qq) in Two Color QCD (TCQCD) and Nambu - Jona-Lasinio (NJL) models (see eg. 
[21 El) inspired by the possibility that the QCD ground state at large \i and small 
T is a color superconductor |3j. However, conceptually simpler issues such as the 
restoration of chiral symmetry and the spectrum of excitations in a dense interacting 
medium continue to present interesting challenges. One obvious example is whether 
and how masses and widths of mesonic bound states such as the p meson are modified 
in a dense medium jH]. A broader question is how a Fermi surface manifests itself in 
Euclidean simulations. Recall that color superconductivity phenomenology is based 
on the BCS mechanism for electronic superconductivity, which in turn is based on a 
small, though non-perturbative, deformation of a sharp Fermi surface. 

As part of the learning process we have studied the simplest non-trivial model 
simulable with |i / using standard algorithms, namely the Gross-Neveu (GN) 
model in 2+1 dimensions. Its Lagrangian density in Euclidean metric is written in 
terms of 4iVf -component spinors ip, ip as 



In the chiral limit m = the model has a global Z 2 symmetry ip i— > 75^, ip 1— > — ^75. 
At tree level the model's content is Nf fermion flavors each with parity-invariant 
mass m, interacting via a four-fermion contact term. However at leading order in 
an expansion in 1/Nf it can be shown that for sufficiently large coupling strength 
g 2 the Z2 symmetry is spontaneously broken by a condensate (ipip) 7^ leading to a 
dynamically generated fermion mass gap E = g 2 (tptp) ^> in. The critical coupling 
g\ at which the gap Hq/Aw —> defines an ultra-violet stable fixed point of the 
renormalisation group at which an interacting continuum limit may be taken. This 
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picture has been verified both at next-to-leading order in 1/Nf and by Monte Carlo 
simulations with finite Nf (eg. [0]). 

The thermodynamic phase structure of the GN model is also known to leading 
order in 1/Nf jZllH]- At T = 0, (tpip) remains constant as /i is increased until the Z 2 
symmetry is restored in a first-order transition at a critical /i c = So. At the same point 
the baryon density n = Nj l {tp^Qip) jumps from zero to nl/ln and then continues to 
rise quadratically as the Fermi disk of relativistic "quark matter" grows. This point 
is hence often called the onset. Simulations with finite Nf have again verified this 
picture with [Sj , and importantly have showed that the first order nature of the chiral 
transition at fi c ^ So persists at small but non-zero T, implying the existence of a 
tricritical point in the (/z, T) plane [H] . With currently available precision there is 
however no evidence for a "nuclear matter" phase with both (ipip) 7^ 0,n 7^ 0. 

The high-/i phase of a related model, the 2+ld NJL model, which has a global 
SU(2)l®SU(2)^ chiral symmetry, has also been extensively studied [10J. The main 
issue of interest has been whether diquark condensation takes place at the Fermi 
surface in the isopsin-symmetric degenerate matter formed for \x > fi c , leading to a 
spontaneous breakdown of baryon number symmetry whose physical manifestation 
is superfluidity. Diquark pairing is favoured in this model because the formulation 
with staggered lattice fermions admits a scalar isoscalar qq operator localised on a 
single lattice site [TT]. It should be noted that in the current Z 2 model a local O(Nf) 
invariant diquark operator is forbidden by the Pauli exclusion principle. Although we 
have not studied the diquark sector in the current work, we will see in Sec. 13. II below 
that the two models show dramatic differences in the spectrum of spin-| excitations 
around the Fermi surface. 

In Sec. I2.1l we will use the technology of the 1/Nf expansion in the so-called "hard 
dense loop" approach to calculate the long-range inter-particle interaction in medium, 
ie. for fi > fi c . The main result is that the Debye mass is infinite, but the plasma 
frequency M a is finite and actually vanishes as the onset \x — > /i c + is approached. The 
main new feature emerging for /i > \x c is the existence of a new physical scale, the 
Fermi momentum hp. Since the lowest energy excitations of the ground state have 
\k\ ~ kp, measurements of Euclidean timeslice correlators with k 7^ are mandatory. 
The resulting energy E(k) is known as the dispersion relation. The main numeri- 
cal results of the paper are thus dispersion relations in various channels of interest. 
In the spin- 1 fermion channel, the low-energy excitations about the Fermi surface, 
which may be hole- or particle-like, are generically known as quasiparticles . There is 
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a well-developed framework for describing quasiparticles and their interactions known 
as Fermi liquid theory. In Sec. 12.21 we will present theoretical predictions of Fermi 
liquid properties to leading non-trivial order in 1/Nf. In Sec. 12. 31 we will describe an- 
other consequence of the new scale; the wavefunction describing inter-particle spatial 
correlations, which decays monotonically in the vacuum, now oscillates with spatial 
frequency related to kp. 

Sec. El describes the results of our numerical study. As mentioned above, the 
fermion dispersion relation is presented in Sec. 13.11 and mesonic (ie. qq) channels 
are similarly investigated in Sees. 13.21 and 13.31 We will find that the predictions of 
the 1/Nf approach are verified with reasonable accuracy, and that in- medium effects 
are clearly visible. In most meson channels the lightest excitations are particle-hole 
pairs with effectively zero energy whose temporal decay is hence algebraic rather than 
exponential, but which nonetheless show a non-trivial dependence on the ratio \k\/fi. 
An interesting phenomenon not described within the 1/Nf approach is a collective 
excitation of the whole system (ie. not a quasiparticle) corresponding to a distortion in 
the shape of the Fermi surface which propagates as zero sound. Tentative evidence for 
an excitation of this form is presented in Sec. 13.31 Finally in Sec. l3.4l we present results 
for wavefunctions in mesonic channels which clearly show the oscillatory behaviour 
described above. After a brief summary, various technical details of the calculation 
of relevant correlators in free-field theory are given in two appendices. 

2 Theoretical Results 

In most theoretical approaches to the GN model it is convenient to encode interactions 
by the introduction of an auxiliary scalar field a related to the fermions via the 
equation of motion a = ■^jipif). Although at tree level a is non-propagating, the 
leading order 1/Nf expansion at fi — predicts that it acquires dynamics through 
quantum corrections due to virtual qq pairs, resulting in a propagator D a (k 2 ) of the 
form [U] 



For soft momenta a thus resmbles an orthodox scalar boson formed as a weakly-bound 
qq state. In the Born approximation a exchange generates an interaction between the 
fermions of range 0(Eq x ) and strength O^Nj 1 ). By contrast, as k — ► oo D a is a harder 
function of k, though still softer than the original tree-level prediction D a oc k°. The 
crossover from IR to UV is encoded as a branch cut in the complex k plane, which 
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means that the decay of D a in Euclidean time is not well-described by a simple pole. 
In the next subsection we will extend the leading order calculation of D a to /x > /x c , 
restricting our attention to soft momenta k <C /x. Even with this restriction the 
resulting dispersion relation is surprisingly complicated. The calculation also gives 
access to the forward scattering amplitude between quasiparticles in the medium, 
which in turn enables contact with the relativistic theory of Fermi liquids ^2] • It is 
hence possible to find quantitative relations between the Fermi liquid parameters kp, 
Bp and /3p, respectively the Fermi momentum, energy and velocity. 

2.1 The a Propagator in Medium 

Here we examine the auxiliary scalar propagator D a in the presence of a non-zero 
baryon density. To leading order in the 1/Nj expansion (eg. |6J, we have 



D-\k) = l-n(*;Ai) 

i r 1 i l 

(3) 



If If 1 1 

tr— h / tr- 



2 

.S Jq,n=0 " H + So ' Jq,n>n c + ^7o Ki ~ JO + ^70. 

where II is the virtual fermion - anti-fermion vacuum polarisation bubble, we have 
used the gap equation at /x = to express 1/g 2 in terms of the zero density gap 
S , and assumed that the gap vanishes in the integral defining II for /x > /x c . It is 
technically more convenient to evaluate the integrals for temperature T > and for 
arbitrary spatial dimension d G (1, 3); after tracing over Dirac indices we find 

d d q 



D-\k ,k) = Ag 2 Tj2 J 



ka(u} n — ko — ifi) + k.(q — k) 



(4) 



[(u n - i/j,) 2 + q 2 ] [(u n -ko- in) 2 + (q- k) 2 ] 

where the sum is over discrete Matsubara modes u n = (2n — 1)ttT. 

The first two terms of (jlj) are individually divergent, but their sum is finite. In 
the limit T — » (see eg. [Hj) their contribution is 

v <f^.f^. y^-fi (5) 



(4,r)fr(|)Uo vWs! A J (*r)Sr(f)(<i 

the critical chemical potential in d + 1 dimensions in the large- Nf limit being 

.HMH,',, (6) 
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The frequency sum for the momentum-dependent term may be evaluated using 
the formulae(eg. [T3j ) 

-tj;avw.- 1 *,«o= £ /"/ff'"^'^ < 7 > 

n s,=±l 

™p »+/ 7-1 \ a — / 771 \ (fJ.-s 1 E 1 )(l-f + (s 1 E 1 )-f-(s 2 E 2 )) 

^— ' 4SiS 2 -C/l-C/2(^ — Sl-G/l — S2-B2) 



Si=±l 



where A ± (u,E) = ((00 =p z/i) 2 + i? 2 ) -1 , the Fermi-Dirac function f±(x) = (exp((x =F 
fi)/T) + and is taken positive throughout. Using the identity 1 — f±(x) — 
/=p(— x) = 0, and noting that for fi > lim T ^ f+(E) = 1 — 9(E — /j l ), lim T ^ f-(E) 
0. we find in this limit 

where the spatial loop momentum has been rerouted by q — > q + 5 k so that £7^2 = 
|g ± We now restrict our attention to ko, \k\ <C /i, thus probing the inter- 

particle interaction for soft momentum transfer in the so-called hard dense loop (HDL) 
approximation 1 . The first and third terms together yield 

V f 2w d0 f°° dqqd -2 k l + \k\ 2 + ik \k\ cos 6 + 0{k 3 cos e/q,k 4 /q 2 ) 



#y~ 3 /^ + |£| 2 .Jfe, 



+ ^|fc| 2 + 0(A;V/i 2 )) (9) 



(47r)fr(f)V 3-d 4// 

The integrand of the second term is non-zero only in a crescent-shaped region near 
the Fermi surface; we find 

2 Ijfel yf d# costf ,„ + iMco S e+... 
<7 - / —77 77; — I ,,, ^99 



/X— ^ cos t 



(47T)lr(f) 7-| 27TA;2 + |fc|2 C OS 2 ^ 
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x (— 2(A;g + |fc| 2 ) + i— |£| 2 cos#sin 2 # + 0(A;7g 2 )) = (10) 



W L (J j + |jE|. )1 + j W 3 , + ^ 3 



(47r)ir(f)*:„ + (fcS + |Sp)H 4^[fe + (fej + V 



2 y 



1 For asymptotically large /; > /i, Sq the form of the propagator reverts to D tJ 1 {k) (x k d 1 
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Eqns. ©, © an d (fTUj) combined yield 



2 d-3 



(47r)lr(f)(3-d) 



4| 3 ^^- d (/- 1 -^- 1 )+ (11) 

*g + W» ( i - (3-^)(4+W) ^ 3 i ) + 

v k + (k 2 Q + \k\ 2 )* ) 4/i v [ fco + + |^|2)i] 2 ; V/i 2 

The first observation to make is that in the static limit ko = the momentum 
dependence of (II lj) vanishes for the physical value d = 2, implying complete screening 
of the static potential due to a exchange, or in other words infinite Debye mass. 
By contrast, for k = D a is proportional to a conventional zero-momentum boson 
propagator with mass 

Ml = 4|^ 3 - V- 1 " ^) (12) 

or for the physical d = 2 

M a = 2 v V(/i-/i c ). (13) 

Just above the transition, therefore, there is a tightly-bound state whose decay at 
rest into a particle - anti-particle pair is Pauli-blocked. The free-field result M a = 2/j 
is only approached as /x — > oo. M CT may also be interpreted as the plasma frequency 
of the system. 

For states in motion, things are more complicated. For |fc| <C \ko\ in d = 2 we 
have 

D~\k , k) cx Ml + k\ + i|fc| 2 + i*^ + 0$jr). (14) 

Eqn. (|14j) yields two pure imaginary poles in the complex ko plane leading to solu- 
tions with differing decay lengths in forward and backward directions in Euclidean 
time. In principle this may be probed by a Euclidean lattice simulation. However, 
linear response theory dictates that on analytic continuation to Minkowski metric, 
the physically relevant solution giving the energy E(k) of a collective excitation of 
the system is given by poles of the retarded propagator , which corresponds to the 
solution of (fT4*|) with negative imaginary part: 

E(i) = {Ml + \W) 1 + f c M. + ^ + -L) . (15) 

In the limit /i — > oo this becomes E ~ M a + \k\ 2 /2M a , appropriate for a free non- 
relativistic particle of mass M a . Another limit we may consider is \i — > /i c+ , in 
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Figure 1: The dispersion relation E(\k\) in units of fi for various values of M a . 



which case <C \k\ 
relativistic boson: 



in (fTTj). yielding the dispersion relation of a standard 



E 2 (k) ~ M 2 + \k\ 



(16) 



Dispersion curves obtained from a numerical solution of (jllj) for three different Mo- 
are shown in Fig. ^ 

2.2 Fermi Liquid Properties 

Given what we know about the fermion - fermion interaction in medium from the 
calculation of D a , it is possible to make quantitative statements in the context of the 
Fermi liquid description of degenerate fermions, first developed for non-relativistic 
systems such as 3 He by Landau |Ej (see also [EH)- The generalisation to relativistic 
systems has subsequently been given by Baym and Chin [T2] . The essential phys- 
ical idea is that the dominant low-energy excitations in the neighbourhood of the 
Fermi surface are long-lived quasiparticles having energy eg, width oc (eg — /i) 2 and 
equilibrium distribution ng related by the Fermi-Dirac function 




(17) 



For T = we expect to have the form 

eg ~ fj, + f3 F 

where k F , (3f = |Vg£|£| =fcF | are respectively the Fermi momentum and Fermi velocity, 
which in the absence of further information are phenomenological parameters. The 
quasiparticle energy is implicitly determined by the equation for its variation under 
a small departure from equilibrium; 

/d d k' 
J^y^P 5 ^ ( 19 ) 

where Ttt, = J-t, t is the Fermi liquid interaction, which is in turn related to the two 
particle forward scattering amplitude 



where we have retained explicit spin states a, a' and KA is the Lorentz invariant 
matrix element. Landau's theory becomes predictive if this amplitude is calculable 
in medium. 




Figure 2: Feynman diagrams for forward scattering, showing direct (left) and ex- 
change (right) contributions. 



To leading non-trivial order in 1/Nf there are both direct and exchange contri- 
butions to Ai, as shown in Fig. but in the chiral limit only the exchange term is 
non-vanishing. After averaging over spin states we find 



k,k> 



1 - 



k.k' 

£ k £ k' 



(21) 
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Specialising to d = 2 and using the result (jllj) for D a we find that for k, k' on the 
Fermi surface separated by angle 9, 

with M 2 = Afi(fi — fi c ). This particularly simple form arises due to the vanishing spa- 
tial momentum dependence of D a for zero energy exchange since the a has infinite 
Debye mass in medium. The factor of Nj 1 ensures that T can be treated systemati- 
cally as a weak correction to the leading order case of non-interacting particles; note 
that this would not hold if direct interactions were involved, since that would result 
in a compensating "degeneracy factor" of Nf in the numerator. 

It is possible to derive relations between the Fermi liquid parameters. For instance, 
by requiring consistency of the above framework under Lorentz boosts the following 
relation may be derived j!2j : 

r d 2 k' 

£ k^k £ k = k + J J^f^k'^k^k'- ( 23 ) 

Dotting this equation with the unit vector k, and noting that at T = 0, n%, = 6>(yU— e^,) 
and Pp5(ep — fi) = S(\k'\ — k F ) we arrive at the following relation at the Fermi surface: 

Hp F + fik F g / — (l-cos#)cos# = k F , (24) 

N f M aW JO 47T 2 

ie. 

k F k F u . . 

"' = 7 +, iW0- (25) 

Here q is the degeneracy of the states participating in the exchange interaction; for 
the Z 2 GN model = 2 counts the number of spin states, whereas = 4 for the 
isospin-symmetric matter described by the NJL model when 7r exchange is also taken 
into account. 

Similarly, the compressibility follows from differentiating the relation \x = Sk F with 
respect to n: 

dn de k dk F - 2ir ir/i 

dn dk F dn + Qk F iV/MJ(//)' 1 ' 

where we have used the equality of particle and quasiparticle densities n = J j0vn^ 
and the bar denotes the average over the circle. Using relation (J23j) we then derive 
the first sound velocity f3\\ 

2 = ndji_ = 3 9 k F 
11 ~ ii dn 2/i 2 8N f M 2 (fi)' { ' 
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To find the relation between e^ F = \i and kp = (Aim/G)? we need to integrate 



F 
kp 



Fc 



1 + 



k F dk F = I - 3g — — . {21 



The integration limits specify the conditions at onset; note that kp c = ^c + Oi^Nj 1 ) = 
S + 0(NJ 1 ) with the corrections given by the two-loop gap equation. The solution 
of (|2~%j) expressed consistently to O^NJ 1 ) is 



16iV 



.f 



M 2 - £ + 2^ - Mc) + 2fi In (3^^) 



(29) 



We can combine this with (j2*3j) to get an absolute prediction for [3 p. For /i ^> ji c we 
find 

- = i-4r; )9f = i--^t, (30) 

which, happily, is consistent with the causality requirement j3p < 1. In fact, the full 
solution (|2T)|) violates this condition for fi — fi c w 0(fi c ), although a precise statement 
is impossible without knowledge of /cp c to 0(NJ 1 ). This suggests that the HDL 
approximation leading to (fTTj) is insufficient at such small \i. 

2.3 Friedel Oscillations 

Another characteristic of a Fermi surface can be exposed by considering spatial cor- 
relations between qq pairs. A convenient way to do this is to consider correlation 
functions of the form 

C(y;x ) = ^(#V>(0,0) ^(x + y, x )) (31) 

X 

(■ (■ p ipx -iqx -iq.y 

= tr / / r- — -r-^ (32) 



p Jq 



where in 2 + 1 dimensions j p = J d 3 p/ (27r) 3 and the matrix T determines the quantum 
numbers of the channel. M is the fermion mass, which we will take to have the leading 
order form M(/i) = £o#(/i c — /•*)• Performing the Dirac trace and the sum over x yields 

4 Li f ,2^(±(Po-^)(go-'i/i)±P 2 + M 2 )e i(p °- 9o)so e"^ 



C{y] xo) = J dq J dp J d 2 p 



[(po - ifi) 2 + V 2 + M 2 ] [(g - %y) 2 + P 2 + M 2 ] 

(33) 

where the ± signs correspond to pseudoscalar and scalar channels respectively. 
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First consider the chirally broken phase p < p c = £ ; t ne Po (Qo) integral may be 
performed using Cauchy's theorem by completing the contour in the upper (lower) 
half-plane, picking up residues at the poles at ip ± iyjp 2 + M 2 . Therefore 

CKO - ^/^-»« x {^_ ,34, 

POO 

= 2 pdpJ (py)e~ 2xo ^ p2+M2 . (35) 
Jo 

where J is the Bessel function and in the second line we have specialised to the pseu- 
doscalar case. For large xq, the final integral over p can now be approximated using 
Laplace's method of asymptotic expansion, by expanding the integrand in powers of 
p/M, yielding 

lim C(y; x ) ~ — e~ 2Mx ° exp ( JJ^M.) . (36) 

The general profile of the wavefunction is a gaussian with width increasing as yjxo. 

In the chirally restored phase p > p c , M vanishes. For \p\ < p,, both poles of 
the integrand of (J33|) lie in the upper half-plane, making the go integral vanish. The 
effect is to raise the lower limit of the integral of (|35|) from to p. The asymptotic 
expansion is now made by expanding the integrand in powers of (- — 1), with the 
resulting form 

lim C(y;x ) ~ ^- e ~ 2 ^ J (^\y\)- (37) 

20-*0 Xq 

The wavefunction profile no longer changes with xq, but instead oscillates with a 
spatial frequency determined by p, which to this order may be identified with the 
Fermi momentum kp. 

The oscillations observed in C(y;xo) are characteristic of a sharp Fermi surface 
and are reminiscent of oscillations in either the density-density correlation function, 
or the screened inter-particle potential, in degenerate systems generically known as 
Friedel Oscillations |16j . The resemblance is only qualitative, however, because while 
the large-?/ form of (|37jl is C oc y~^- sia(kjry— |), the form of eg. the screened potential 
in high-density QED is y~ 3 cos(2kp y), and even y~ 4 sin(2kpy) in the relativistic limit 
n^e *C p [El • The disparity is due in part to the lower dimensionality, and partly be- 
cause the function C(y; xq) probes correlations between fermion field variables rather 
than scalar densities. In Appendix A we present a similar asymptotic analysis of 
the density-density correlator C nn (y;x ) = ^0' (0, 0)jj(x, x o )j {0, 0)jj(f + y, x )), 
where jo(x) is the local baryon density ^^^(x), and find a decay with y more closely 
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resembling the forms in the literature. In particular, because C nn describes cor- 
relations between particle-hole pairs at the Fermi surface, the asymptotic spatial 
frequency of the oscillations is 2kp rather than hp. Because of the similar origin, 
however, we will refer to oscillations of the form (J3*7j) as Friedel oscillations in what 
follows. 

3 Numerical Results 

In this section we present results from a numerical simulation of the Z 2 GN model 
using staggered lattice fermions. The lattice formulation and simulation algorithm 
are described in detail in Refs. [HUE]- We used a 32 2 x 48 lattice and throughout chose 
a bare fermion mass m = and N = 2 lattice flavors. Unless otherwise stated the 
coupling constant used was a/g 2 = 0.75, corresponding to a physical fermion mass at 
/i = of S a ~ 0.17. At each value of /i between 2000 and 4000 hybrid Monte Carlo 
trajectories of mean length 1.0 were taken. All results have been taken in the chirally- 
restored phase with /i > /i c , whose value is estimated to be \i c a ~ 0.16(1). Values for 
the baryon density per staggered fermion flavor na 2 for values of chemical potential 
/ia G [0.2, 0.8] are given in Table [TJ showing that the density rises steadily as \i is 
increased. The reported densities are roughly twice the free-field value n = [i 2 /2ti, 
consistent with fermion doubling; N flavors of staggered lattice fermion describe 
Nf = 2N continuum flavors in 2+1 dimensions, so that Nf = 4 in our simulation. 
Henceforth we will quote results in units in which the lattice spacing a — 1. 

Table 1: Baryon density n = N~ l (i/j^qi/j) 





n 


0.2 


0.0113(1) 


0.4 


0.0520(1) 


0.5 


0.0865(1) 


0.6 


0.1394(1) 


0.8 


0.3127(1) 



3.1 Fermion Dispersion Relation 

The fermion timeslice propagator at fixed spatial momentum k is defined by Q (k, Xq) = 
J2x&(®>fy%' x o) e ~ lkx > where Q(x,y) = (ip(x)ip(y)). The decay in Euclidean time xq 
then yields the fermion dispersion relation E(k). For a degenerate system with Fermi 
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momentum kp, excitations with \k\ < kp are hole-like and those with \k\ > kp 
particle-like; the lowest energy excitations are those in the neighbourhood of the 
Fermi surface with \k\ ~ kp. Fig. |3]shows E(k) measured for k parallel to the x-axis 
for four values of \i > fj, c , using the procedure specified in [TU]. In order to obtain a 
smooth curve we plot as negative for \k\ < kp. Note also that \k\ is plotted as a 

fraction of 7r - since the staggered fermion action is only invariant under translations 
of length 2a the maximum momentum which can be probed is it/2. 




Figure 3: The fermion dispersion relation -E(|fc|) at four values of \x. 

It is obvious that there is no mass gap at the Fermi surface characteristic of a 
BCS instability, since for certain combinations of fi and k our data have ~ 0. 

Indeed, it is possible to fit a smooth function 

E(\k\) = -E + D sinh-^sin \k\) (38) 

to the data, which corrects the Fermi liquid form (J18)) for lattice discretisation effects, 
and permits the extraction of an effective Fermi momentum Kp and velocity (5p given 

by ma 

^ = sinh- 1 (sinM = ^; (3f = D^^. (39) 

D cosh K F 

The fitted values of Kp,/3p are given in Table 121 The departures from free- 
field values are small and of similar magnitude to the predictions (}30|) of the 1/Nf 
expansion Kpj\i ~ 0.953, /3p ~ 0.984. Also shown is the ratio Kp/[i/3p, which apart 
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Table 2: The Fermi momentum and Fermi velocity. 






K F 


Pf 


K f /hPf 


0.2 


0.18(1) 


1.04(2) 


0.87(5) 


0.4 


0.386(6) 


0.96(1) 


1.00(2) 


0.6 


0.584(3) 


0.98(1) 


0.99(1) 


0.8 


0.782(14) 


0.96(1) 


1.02(2) 



from jj = 0.2 shows no significant deviation from its free- field value 1. This constrasts 
sharply with Monte Carlo studies of the 2+ld NJL model, for which q — 4, where 
values in the range 1.3 - 1.5 were observed [TU]. This suggests that the quasiparticle 
spectrum observed in that model is determined by physics outside the scope of the 
I/TV/ expansion, and as such the current result provides indirect support for the exotic 
gapless superfluid scenario proposed in [TU] . 

3.2 The a Dispersion Relation 

Just as in the continuum 1/Nf expansion, the lattice formulation of the GN model 
encodes the interactions between the fermions via an auxiliary bosonic scalar field a 
whose value is given by the equation of motion o = -j^-^^ (in this paper we define 
cr with dimension |). For /i = and fc < E the auxiliary propagator (j2J) resembles 
that of a standard boson of mass 2So. Since this mass coincides with the threshold 
for an qq continuum, D a is dominated by a branch cut rather than an isolated pole 
in the complex fc-plane making standard lattice spectroscopy techniques difficult to 
apply A recent study using the Maximum Entropy Method has shown evidence 
for a non-zero binding energy in the a channel for finite Nf [T5] . 

For n > fi c , however, the calculation of Sec. 12.11 suggests that the a is more 
tightly bound. We have studied this by measuring the auxiliary timeslice propagator 
V(k,\x \) = Yl,3 a;o=±|a;o| 0)cx(x, x ))e~ , once again choosing k along the x- 
axis 2 . Fig. HI plots V(xo) for k = 0, showing that with the exception of ji = 0.2, 
simple pole fits of the form T>{xq) oc e _M<TX ° describe the data reasonably well before 
it descends into noise. Values of M a extracted from fits to timeslices Xq G [1, 5] are 
plotted in Fig. EJ Also shown is the leading order prediction M a = 2^J ^{^i — /i c ), 
with a value 0.16(1) assumed for /i c . The numerical results fall 20-30% lower than 
the leading order prediction, which is consistent with a correction of Oi^Nj 1 ). 

Results from simple pole fits to T>(k, \xo\) for \k\ < | yielding the a dispersion 

2 Note that as discussed below eqn. (|14f) T> may not in fact be time-reversal symmetric for k ^ 0. 
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Figure 4: The a timeslice correlator T>(xq) at k = for 4 values of \i. 




jure 5: M a as a function of fi. The line is the leading order 1/Nf prediction 
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Figure 6: The a dispersion relation E{\k\) for 4 values of fi. 

relation E(k) are plotted in Fig. El 3 . Only for the lowest value \i = 0.2 is there any 
marked change in E as k is increased, in which case E ~ k. For larger \x the dispersion 
is flatter. There is qualitative agreement with the leading order 1/Nf result shown in 
Fig. HI although an improved understanding of discretisation effects close to \k\ = | 
would be needed to make the comparison more quantitative. Nonetheless, Figs. |U 
and |H1 demonstrate unambiguously that in- medium effects are observed in the lattice 
simulation. 

3.3 Mesonic Dispersion Relations and Zero Sound 

Next we investigate mesonic (ie. qq) states. We have studied mesonic correlation 
functions Cr(k, Xq) = X^^0'r(0, 0)jp(x, xo))e~ lk,x where the bilinears jr(^) are defined 
with scalar, pseudoscalar or vector quantum numbers; in terms of staggered fermion 
fields x, X the operators are 

ji{x) = XxXx] j 15 (x) = e x XxXx] hi(x) = -^-[XxXx+i + Xx+iXx], (40) 

3 Since under the global Z2 <j(x) (— l) Xo+Xl+X2 <r(x), in the symmetric phase the propagator 
obeys V{ki) = V(tt — fc,). 
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(1=0.2 



(1=0.4 




— * — * 

Figure 7: Pseudoscalar correlator C 75 (\k\,Xo) at 4 different fx for momenta \k\ = 
(filled circles), ^ (filled squares), | (filled diamonds), || (filled up triangles), j (filled 
down triangles), jj? (open circles), ^ (open squares), || (open diamonds) and | (open 
up triangles). 

where = (— l) x °, r) 2x = (—l) Xo+Xl and e x = (— \Y 0+Xl+X2 - j± s De f ore we study the 
timeslice correlators in each channel as a function of spatial momentum k \\ x. 

As expected for ii > fi c where the Z 2 global symmetry is restored, the results in 
scalar and pseudoscalar channels coincide. The results on odd timeslices, averaged 
over k and it — k, are shown for four different \i > fi c in Fig. [7| Note the difference 
in the vertical logarithmic scales between the upper two panels /i = 0.2, 0.4 and the 
lower two fi = 0.6,0.8. For each /i, C 75 (/c = 0) is almost flat, indicating a massless 
state. 

— * 

Appendix B sketches the calculation of Cr(k,xo) as a function of /i in free field 
theory. As Xq — > oo the function is dominated by a continuum of particle-hole pairs 
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at or near the Fermi surface, which effectively cost zero energy to excite. The generic 
result is that for \k\ < 2/i the decay is algebraic, with ()60J67|) 



j x i 


Ijfel 






[ X 2 5 


\k\ 



C r (k,x )oc{ Ta' Z]^^ (41) 

2/i. 

Only once \k\ > 2/i does it become kinematically impossible to excite a pair with zero 
energy, resulting in exponential decay: 

C r {k,x ) oc x 2 exp(-(\k\ - 2p)x \ (42) 

The sequence of plots in Fig. [7| is in qualitative agreement with these findings. For 
each /i there is a particular value of \k\, highlighted with a solid line in the plots, 
for which the temporal falloff is particularly slow, corresponding to \k\ ~ 2/i (eg. 
for /i = 0.4 the slow falloff occurs for \k\ — ^ ~ 0.785). For \k\ larger than this 
value the decay is much steeper, although only for /i = 0.2 does it resemble an 
exponential form. Because of the technical difficulties in treating correlators with 
power-law decays in a finite volume (see eg. 0) we have made no attempt to fit the 
numerical data for C 75 (|fc|, xq) to a functional form. It can be observed, however, that 
the overall magnitude of C 75 increases with p as the size of the Fermi surface and 
hence the number of participating particle-hole states grows. There is, however, no 
obvious correspondence with the prefactor \k\/p predicted in (JHUJ), probably because 
the approximation \k\ <C /i is not valid over the accessible momentum range. 

The situation is more interesting in the vector channel, corresponding to the quan- 
tum numbers of the p meson, because for \k\ > we can distinguish between C 7| (&), 
in which the component of the vector is parallel to k, and C JX (k). In Fig. |H]we plot 
the correlator for several \k\ values at p — 0.6 in each case. For C 1± , just as for 
C 75 , there is no difference between the data for \k\ and it — \k\. For C 7|| , however, 
because of the point-split nature of the operator j 7i this symmetry no longer holds; 
in this case we choose to plot data with \k\ G [§,7r] (the data for \k\ < ~ are very 
similar to those for C 7± ). Each plot covers two decades on the vertical axis; however 
both the magnitude and the shape of the curves is very different. For C 7x the curves 
are qualitatively very similar to those of Fig. at p — 0.6, with a distinguished 
momentum \k\ = The correlator C 7| is much smaller in magnitude, and is con- 
sistent with exponential rather than algebraic decay. The lines are fits of the form 
C 7|| (|£|,x ) = A{e~ Exo + e -^{Lt-x Q )y The resulting E ^ is shown in Fig _ ^ For 

small \k\ it resembles that of a massless pole, ie. E = Po\k\, with velocity (3q ~ 0.5. It 
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Figure 8: Vector correlators C 7± (\k\, xo) (left) and CyJir — |&|,£o) (right) at p = 0.6. 
The symbols have the same meaning as in Fig. 

should be mentioned, however, that although C 7|) and C 7± still differed, no evidence 
for a massless pole was seen in the data at p = 0.8. 

Light states in the p channel in medium have long been of interest [3] . One recently 
proposed theoretical scenario is that the p should be viewed as an "almost" Goldstone 
boson due to a spontaneous "induced" breaking of Lorentz symmetry over and above 
the explicit breaking due to p ^ jTH] . In this case the p mass is predicted to scale as 
m 2 p oc 2pn. Another, proposed in the context of Two Color QCD (TCQCD) in which 
mesons and baryons fall in the same multiplets and are degenerate at p = 0, is that 
for sufficiently large p there is a condensation of vector diquark states leading to the 
spontaneous breaking of rotational invariance and vector Goldstone states with both 
E oc \k\ and E oc \k\ 2 j20J - In- medium modification of the p correlator has recently 
been observed in lattice Monte Carlo simulations of TCQCD |21j. 

Since our results suggest a massless state and there is no Fermi surface in TC- 
QCD we favour an alternative explanation, namely that the pole signals a collective 
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Figure 9: Dispersion relation E{\k\) extracted from C 7|| . 

excitation of the system in which the shape of the Fermi surface is distorted, leading 
to a rotationally non-invariant disturbance propagating at velocity (3$ known as zero 
sound. The conditions for zero sound to propagate are that there must be interac- 
tions, ie. J-^^i 7^ 0, and that inter-particle collisons are negligible, ensured if T ps 
|22[I15|. Indeed, zero sound solutions can in principle be found self-consistently using 
the theoretical framework of Sec. l2.2| although the specific form (|2*2*|) for the Fermi liq- 
uid interaction entails the solution of an integral equation of the second kind, beyond 
the scope of this paper. A potential problem for this picture is that the extracted 
velocity (3q ~ 0.5 is less than (3p ~ 1, implying that quasiparticles would experience 
damping via Cerenkov radiation of zero sound. Clearly further work exploring the 
systematic effects of varying /i, Eo, k and volume will be needed for a more complete 
understanding to emerge. 
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Figure 10: Scalar wavefunction at /j = 0.5 for (5 = 0.65, 0.70, 0.75. 

3.4 Mesonic Wavefunctions and Friedel Oscillations 

We now switch attention to spatial correlations between particles, probed via the 
wavefunction ty(x) defined by 

* r (x) = lim ^\x = 0) J2(G q (0, 0; y, x )Tg q {0, 0;y + x, x )T) (43) 

y 

where as ususal T projects out the quantum numbers of the channel of interest. In 
the GN and related models ^>(x) is technically much easier to define and measure 
than in Q CD-like theories where they are gauge-dependent. Meson wavefunctions in 
the 2+ld GN model have been studied at T, fi — in where further technical 
details are given. 

Fig. Hill shows ty(x) measured along the x-axis in the scalar channel at a constant 
ji = 0.5 at various values of the the coupling 1/g 2 = 0.65,0.70,0.75. The oscillatory 
behaviour described in Sec. 12.31 is clearly seen, and is plainly not a discretisation 
artifact since its form is stable as a(g) varies. Fig. II II shows that there are no signifi- 
cant differences among the various mesonic channels implying that in contrast to the 
situation at /x = [2H], effects due to eg. a exchange are very hard to detect. 

In the next four figures we plot the scalar wavefunctions at 1/g 2 = 0.75 for a 
sequence of /i values 0.2, 0.4, 0.6, and 0.8. The simulation data are connected by 
solid lines. As /i, and hence Icf, increases the oscillations decrease in wavelength 
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Figure 11: for scalar, pseudoscalar and vector channels at 1/g 2 = 0.75, \i = 0.5. 



in accordance with theoretical expectation, and provide a graphic illustration of the 
presence of a sharp Fermi surface. The dashed lines show measurements taken with 
the same lattice parameters but with the interaction switched off. The disparity with 
the interacting theory is small, though increasing with |x|, showing that the free field 
description of the oscillations is qualitatively correct. A possible explanation is the 
infinite value predicted for the Debye mass in Sec. 12.11 implying that interactions 
between static quarks at non-zero \x\ are completely screened. Also shown is the 
theoretical form (|S7|) = J (k F x) with k F = /i, showing good agreement with the 
data for small x ^5 2kp 1 . Unfortunately it appears hard to obtain more quantitative 
information, such as an independent fit for the Fermi momentum kp, because \Jo(kx)\ 
decays only as x~z. This means that fits should not only include the backwards- 
propagating signal J (k F (L s — x)) but also image contributions J^kpinLa — x)) [S] - 
our attempts to find a satisfactory fit were unsuccessful. The figures therefore simply 
show both the "forwards" Jo(fix) and "forwards-and-backwards" Jo(fix) + Jo(fi(L s — 
x)) forms, showing that neither gives a satisfactory description of the data over the 
full range of x and /i. It is also possible that more sophisticated fits taking proper 
account of discretisation effects are needed; at \i = 0.8 the baryon density n = 0.31 
(see Table [Q, a significant fraction of its saturation value of 1, so that the Fermi 
surface may be distorted. 
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Figure 12: Scalar wavefunction at ji — 0.2. 




Figure 13: Scalar wavefunction at ji — 0.4. 
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Figure 14: Scalar wavefunction at ji — 0.6. 




Figure 15: Scalar wavefunction at ji — 0.8. 
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4 Summary 



The GN and related models remain the only interacting field theories both simulable 
by standard lattice methods at fi ^ and displaying a Fermi surface, thought to be an 
essential feature of dense quark matter. In this paper we have attempted to develop 
an understanding of how orthodox lattice observables are altered in such conditions, 
and how they should be interpreted using the language of many-body physics. The 
main achievements have been: 

• An analytic calculation of the auxiliary boson propagator D a in medium to lead- 
ing order in Nj 1 in the HDL approximation. The branch cut in the complex- 
plane is modified to become an isolated pole. Physically, the Debye mass is infi- 
nite, but the plasma frequency finite and vanishing as onset is approached. The 
dispersion relation E(k) is non-trivial. Lattice simulation has verified these re- 
sults with acceptable precision (ie. to within 0(NJ 1 )), and have unambiguously 
identified in-medium effects in this channel. 

• A systematic calculation of Fermi liquid parameters up to O^Nj 1 ), which is 
consistent with causality constraints and has in turn been verified by lattice 
simulation. This contrasts with existing studies of the quasiparticle dispersion 
relation in the 2+ld NJL model with \x > fi c , which are not consistent with the 
large- Nf predictions. 

• A study of meson correlation functions at k ^ 0; instead of showing exponential 
falloff with Euclidean time these generically decay algebraically, signalling the 
presence of massless particle- hole excitations for all k. There is qualitative 
agreement between the results of free field theory and lattice simulation, which 
both show a non-trivial dependence on \k\/ fx. There is also tentative evidence 
for a massless pole in the vector channel, which is possibly a manifestation of 
zero sound. 

• The first observation of oscillatory behaviour in mesonic wavefunctions, which 
resemble the Friedel oscillations familiar in many-body physics. Figs. ll2flT51 offer 
a graphic confirmation, if one is still needed, of the presence of a sharp Fermi 
surface in this model. 

It would be interesting to apply the same analysis to the lattice NJL model at 
ji 7^ 0, which has been studied in both 2+ld [H] and 3+ld |3j- Independent of 
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the issue of whether the ground state has a non- vanishing diquark condensate, these 
models are more realistic because they contain physical pions, whose behaviour at 
the quark-hadron transition is thought to change dramatically. Indeed, the result 
of Sec. 12.21 that the plasma frequency (pole mass) in the chirally symmetric phase 
decreases while the Debye (screening) mass diverges is reminiscent of 

recent results obtained by Son and Stephanov in the chirally broken phase j21] using 
an effective theory and scaling arguments. It would be straightforward to repeat 
our analysis for the pion dispersion relation. More generally, it is worth noting that 
should an algorithm effective at fi ^ ever emerge our results in the meson sector 
are in principle accessible as gauge invariant correlators in a QCD simulation. The 
a channel in the current model is dominated by disconnected qq bubbles which are 
technically hard to compute if an auxiliary field is absent, but in the lightest channel 
in QCD, namely the pion, the relevant diagrams have connected quark lines. The 
main difference in the baryon sector, of course, is that gauge invariant quasiparticles 
now have qqq quantum numbers, and at least below the transition should be non- 
relativistic at the Fermi surface. 



Acknowledgements 

S JH was supported by a PPARC Senior Research Fellowship, CGS by the Leverhulme 
Trust and DOE grant DE-FG02-96ER40945, and JBK by NSF grant NSF-PHY- 
0102409. We are grateful to Kurt Langfeld for discussions in the early stages of this 
work. 



Appendix A 

In this appendix we outline how oscillatory behaviour as a function of spatial separa- 
tion y resembling Friedel oscillations appears in the large Euclidean time xq — > oo limit 
of the density-density correlator C nn (y; x ) = X^O'o(0, 0)j j (x, x )j (0, 0)j j (x+y, x )), 
where jo(x) = ip r Y ip(x). The relevant expression is 

f T f f \ e ipX e ikx e ik.y ^-iix e ~ii.y | 

2Z tr / / / M^TT To-:— 7ott— 70-7- >, (44) 

~ J P J q Jk h\ % V + ^70 W + /^7o if + At7o 4 + /^7o 
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where we have assumed massless quark propagation in the high density phase. The 
trace yields three terms of which the first is 

\p4~ (Po - ~ in)}e i{p - q)x f f [k.£- (k - ifi)(£ - i^e^-i)*^-?)* 



p J q 



[(po - %lif + p*][(q - + q 2 ] J k J e [(k - + k*\ [{io-W + £ 2 } 

(45) 

The integrals over timelike momenta are performed by completing the contour in 
the upper half-plane for f dp , f dk and in the lower half-plane for f dq , J d£ , 
noting the pole positions at eg. po = in ± i\p\. In fact, only the pole with negative 
imaginary part contributes to J dq and J d£ , and it is straightforward to show that 
the dominant term for large xq in the resulting expression comes from the poles at 
Po = i\i — i\p\ and k = ifi — i\k\ since all others are suppressed by powers of e -Atx °. 
The resulting contribution from all three trace terms is 



1 f d d p [°° d d q r d d k [°° d 



e i(p-q+k-e).x e -(\q\-\p\ + \l\-\k\)x 0e i(k-i).y x 



(2tt)^ (2ttYJ (2tt)^ {2v) d 
{(1 + cos8 pq )(l + cosQm) - (1 - cos# pfe )(l - cos6 qi ) + (1 + cos6^)(l + cos 6^)} ,(46) 

where 6 pq denotes the angle between vectors p and q, etc. 

Now, the sum over the timeslice Ylx resu lts in an overall momentum-conserving 
delta-function. Since, however, the integrand is dominated by momenta in the im- 
mediate vicinity of the Fermi surface, a physically reasonable approximation, which 
preserves all angles under the remaining J d d p J d d q, is given by 

(27i) d 5 d (p - q + k- I) ~ {2ir) d 5 d {p + k)(27r) d 5 d {q + £). (47) 

The expression fj4^j) reduces to 

f (g.^-l^-eW^^-D. (48, 

Denote the angle between p (q) and y by d\ (9 2 ). Then 

cos 9 pq = cos 6\ cos 62 + sin Q\ sin 62 cos Q pq , (49) 



where Q pq is the angle between the vectors projected onto the subspace Cl orthog- 
onal to y. The angular integrals over this subspace may be performed using 

d Q<i-i = Zn ■ dQ d - 1 cos 2 Q= , n , (50) 

r(^fi) '7 d-i r(^=i) ' 1 ; 
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to yield 

i 2tt^ 1 r 



(27r) M r 2 (^i) Jo 



PfJ, POO P7T /*7T 

/ p^dp q d - l dq e- 2{q ~ p)xo d6 1 d9 2 e- ipycosei e iqycose2 x 

sin**" 2 0! sin^ 2 # 2 jcos 2 6 1 cos 2 # 2 + -^-j sin 2 X sin 2 # 2 - 1 J .(51) 



The remaining angular integrals are readily performed to yield 

^ ~ ]] p d ~ l dp f q^dq e-^ q - p)xo x 

2(2tt) Jo J^ 

— ] <^ J*(py)J±-i(qy) J±(qy)J^(py) H ^Jd{py)Jd(qy) [(52) 

J i py 2 2 91/ 2 2 2 2 



The integrals over vector lengths p and g are estimated using Laplace's method; 

rli poo ij2d—2 



hm / dpp d ~ l pe 2pX °f{py) / dqq d -\-^g{qy) ~ "f^gfry) (l + O(x x )) , 

(53) 

to yield the final result 



2 d 



(d-l) /i 2d_2 



One may now use the asymptotic form of the Bessel function to deduce the long- 
distance behaviour of C nn (\y\;xo): 



1 cos 2[iy 



, d = 2; 



Cnn(y;Zo) ~ <{ 16 ^ g0 Bin y 2 W ' ,_o ( 55 ) 



Appendix B 



In this section we analyse the behaviour of mesonic correlators as a function of spatial 
momentum Euclidean time in the high density phase. For simplicity we fix d = 2. 
The general expression is 

r (h- r \ = V" ! ! tr ( r HP-7 ~ i(Po - ifJ.ho)T(-ig.7 - i(g - ¥ho)} i(p - q ) x 
r[ ' o) ^]J q [(Po-^) 2 + b1 2 ][(go-^) 2 + |g1 2 ] ' ' 

(56) 

where we will consider Y e {1, 75, 7||, 7±, 7 }, || and _L denoting components parallel 

— » 

and perpendicular to k. On performing the trace, the timeslice sum and the 
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timelike momentum integrals, we are left with, eg. 



C 15 (k;x ) 



e^-\p\)e(\p + k\-^) - - — - -e 

(2*0 \p\\p + k\ 



(57) 



where we have assumed Xq > and once again focussed on the contribution from the 
dominant pole po = i(fi — \p\). We can evaluate the remaining integral in two limits: 

\k\ <C fi: 

In this case the integration region is a narrow crescent in the vicinity of the Fermi 
surface, and the integrand and limits can be expanded in powers of |A?|/|p|- The 
general form is then 



f fiQ 

C r (k;x ) = 2fi\k\ I —h r (9) exp(-\k\x cos9)(l + 0( — 



* 2vr 



where the leading contribution for each channel is given by 
h 



"7s 



\k\ 2 

-hj = - — -cos^sin 2 ^; h. 



4/i 2 



7± — cos 9 sin 2 9; h JQ 



cos#; 



(58) 



h. 



3\k 



k\ 



cos 6 9 + -p- 1 cos 2 9 sin 2 9 + y- cos 9 sin 2 9(1 + 2 sin 



(59) 



■ 2// ■ 4/i 2 

The angular integral leads to the following results for the leading behaviour in \k\/fJ,: 



C 1Bt i(k; xq) 
C 1± (k;x ) 
C 7| | (k; Xq) 



27T/1 



■ ri + n (L 2 (\k\x ) - I- 2 (\k\x ) 



C 75 (k;x ) 



\k\x 
2j! +0(x-) 



\k\ 



2h/j,Xq 



0(xq 



Ti\k\x1 



l*l J 

ll|fc| 3 | 3\k\ 6 s L 3 (\k\x ) - I- 3 (\k\x ) 
307r/i 2/i V (|jfc|a;o) 2 
u \k\ 3d \ (L 2 (f3) -I- 2 (P)\ d 2 



(60) 
(61) 
(62) 



+ lkl [ty 2dp 







3\k\ 



^ / W i FTi _4\— 1 /|D3„5\-1 



C 7o (A;; 



2vr/ia;g vt^x 3 , vr|fc| 3 x 
2//|jfe| 



O (/x|fc|x 



/3=|fc|x 

5J\-1 



7T 



l + -(L 1 (|A ; |xo)-/-i(|A ; |xo; 



2 " +0(x-) 



7T | | OTq 



(63) 



where I and L are respectively modified Bessel and Struve functions. Fortunately 
their difference is a tabulated function [22], which enables the asymptotic form of 
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the second equality in eqns. ([BUI - lo^j) to be determined. The decay in the timelike 
direction is generically oc Xq 2 , the interesting exception being for T = 711 where there 
are significant corrections at shorter distances. 

\k\ > 2fi : 

In this case the 9 functions in (joT|) yield unity for all \p\ < fi, simplifying the 
integration limits. We will restrict analysis to the case T = 75: 



2~ 



fj_ 

71 



7T [ P d P 
^ Jo 



du I dvv 



p 2 + p\k\ cos 9 



-px 



1 i 2fc Q 

1 + — cos 9 

V 



1*1 2 



p 2 (l + -^COS^ + Mr) 2 
1 + 22±* u 2 

W — K 

" [(l+« 2 )(l + 2 « 2 )]^ 



(i>-re) i + (t)+K) 2 « i 



1+ti 2 



— I' 



1 + M 2 



(64) 



where in the second line we have used the substitutions u = cot \9, v — p/fi, and 
k — \k\/ jj, > 2. To focus on the large-Xo behaviour we now employ Laplace's method, 
expanding the integrand in powers of u 2 and noting that since v — k < the negative 
branch of the square roots should be chosen. The dominant contribution is 



dvve -HX {K-2v) 



00 „ ,,2 

2^ u 2 /i 



due 



-flXQ- 



1-) 



•«^o / dvv^{K-v)h 2>1X0V . 
Jo 

(65) 

The final integral is estimated by replacing the factor (k — v) 2 , which is non-singular 
within the integration range, by (k — 1)2; we obtain the lower bound 



C 75 («; s ) > 



2(k-1) 



(66) 



where 7 is the incomplete Gamma function, whose asymptotic form for large fiXo 
yields the final answer 



C l5 (k;x ) > 



^2 e -(\k\-2^)x 



(l + O^xo)" 1 



(67) 



For \k\ > 2fi the correlator thus decays exponentially; however for the critical value 
\k\ = 2 fx the decay is algebraic, and moreover slower than that for \k\ <C fi (JfJDJ). 
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